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DYNAMIC STALL: AN EXAMPLE OF STRONG INTERACTION BETWEEN 
VISCOUS AND INVISCID FLOWS by Jean- Jacques Philippe 

Office National d'Etudes et de Recherches Aero spa t ia 1 es (ONERA) 
(National Office for Aerospace Studies and Research) 

92320, Chatillon, France 

This paper reports on the phenomena associated with airfoils in a dynamic 
stall configuration with emphasis on those associated with airfoils 
oscillating in pitch. The most characteristic experimental results 
of flow separations of '-a vortex nature and their impact on local pressures 
and total' forces are analyzed. A few aspects Tof themethods for pre dict ing 
flows in the presence (or absence) of boundary-layer separation are 
examined, as well .as the main simplified methods available to date for 
computing the total forces in such configurations. 


1 - INTRODUCTION 

The term "dynamic stall" includes a group of aerodynamic phenomena which 
take place when an airfoil', a’ blade section, a- wing or an aircraft is 
subjected to time-varying aerodynamic conditions which lead to drops in 
lift or .sudden increases in the pitch moment which are two characteristics 
of stall configurations. 

For example, during certain helicopter flying configurations, the rotor 
blades can be periodically subjected^to dynamic 'reattachment and stall. 
Indeed, incidence angles on the receding rotor blade can be far greater 
than the airfoil steady stall angles whereas incidence angles of •- the * 
forward-moving blade are low. 

To illustrate the major-phenomena associated with a change of incidence 
angle of the- airfoil about its stall incidence angle, Figure 1 extracted 
from reference 1 shows that,' by comparison with steady aerodynamic 


- 1 - 



tharacteristics , the airfoil can: 

- reach lift levels much higher than those obtained in a stationary flow 
while being the center of major separations with organized vortices shown 
in the figure by means of smoke released at various points on the airfoil, 

- display, during the critical phase when the airfoil is subjected to 
separations, very high-level pitching moments (moment stall preceding lift 
stall) , 

— be subjected to aerodynamic forces sujeh that the work by these forces 
having become positive tends to make the incidence motion unstable. 

In what follows, we will analyze the various aspects of this phenomenon 
and we will cover what we know on: 

- unsteady viscous flows with separation and reattachment over an airfoil 
moving cyclically, 

- the consequences of -such flows on loca;l pressures and on total unsteady 
aerodynamic forces, 

- computation methods for flows about the airfoils that can be applied 
to dynamic - stall' configurations, 

- -practical methods that manufacturers can use to evaluate aerodynamic 
performances under such conditions. 

The author was greatly helped by W.J. McCroskey and his.'-USAAMRDL team 
from the Ames laboratory in the preparation of this document, thanks to 
the work performed jointly in ,the area of dynamic stall and to summary 
reports on this subject (references 1 and 2). Let them be thanked here, 
for their assistance. 

2 - ANALYSIS OF THE PRIMARY CHARACTERISTIC OF DYNAMIC STALL 
2.1. - The vortex nature of separation ijlows. 

Displays of flows about airfoils make itj possible to analyze the type and 
structure of the flows on which depend the aerodynamic characteristics 
of the airfoils. ONERA has a hydrodynamic tunnel which proved to be a 
very useful investigationJ:ool"for studies in unsteady aerodynamics 
(reference 3 ) . • 
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PITCHING MOMENT NORMAL FORCE 



0 5 10 ' 15 20 25 

INCIDENCE, a, deg 


Figure 1 - Example of the development of forces, and pitching moments 
in steady and unsteady'- flow. 
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Figure 2 shows the NACA 0012 profile with pitch harmonic oscillations 
on both sides of the incidence angle for .steady flow stall. When the 
incidence angle increases, we can greatly exceed this incidence without 
seeing the separation that is noticeable on the airfoil (photo 2a). 
However, there is a time in the cycle when large separations take place; 
these display a highly vortex nature. As we can see in Figure 2b, a 
vortex formed near the leading edge moves toward the trailing edge of 
the airfoil. 
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Figure 2 - Pitch Harmonic Oscillations. 

• a = 12.5° + 5° sin'wt; k=0.95; Re = ,6000 

, c 

Displays obtained at the ONERA .hydrodynamic tunnel 
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figure J shows the same axrtoxl set to a txxed geometric xncxdence angle 
close to the unsteady-flow stall angle but wi'th harmonic oscillation beats 
(vertical translations of the airfoil). The rising phase of the airfoil 
retards the appearance of stall on the airfoil as during the phase of 
increasing incidence of the airfoil oscillating in pitch (photo 3a). 

When separations appear on the airfoil, they also have a vortex character. 
We then observe, from a point close to the leading edge, the formation 
and then the upstream motion of' a' v vortexv(Photograph 3b). 
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Figure 3 - Beat Harmonic Oscillations (Up-and-Down Motion) 

<x =12.5°; Az - 0.25 c sin at: k = 0.95; Re ~ 6000 

. c 

Displays obtained at the 0NERA hydrodynamic tunnel. 
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Figure 4 still involves .the same profile set to a fixed geometric incidence^ 
•angle but it is driven through a harmonic back and forth motion (horizontal^ 
translations, of the airfoil). The displays .show that the separated area 
over the top of the advancing airfoil is the center of a continuous source 
of vortices carried away by the flow (Photograph 4a) while, in the case of 
the receding airfoil (Photograph 4b)', the vortex nature is more diffuse. 
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Figure 4 - Back-and-forth Harmonic Oscillations. 

a = 12.5°; Ax = 1.5 c sin ut,k = 0.1075; Re ~ 6000 
5 - c 

Displays obtained' at the ONERA hydrodynamic tunnel 


These displays made in water involve low Reynolds' numbers' (at the ‘chord) 

• i 

in the order of 6000, ■ The USAAMRDL of the Ames Laboratory has a hydro- 
dynamic tunnel that makes studies possible with Reynolds numbers in the 
order of 0.5 x 10^ (Reference 4). Displays are made using hydrogen hubbies 
produced through electrolysis using a technique developed by K.W. McAlister 
of the US Army AMRDL of Ames. The photographs in Figure 5 extracted from 
a document to be published (Reference 5) show the good quality of the 
displays obtained in that manner with a modified NACA 0012 airfoil 
oscillating in pitch, - 
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, 3 


Figure 5 - Displays using hydrogen bubbles (Ames USAAMRDL) . 

NACA 0012 airfoil oscillating in pitch, 
a - 15° + 10° sin a) t; photographs taken with a = a Max 


Displays made in the air by generating smoke from different points on 
the airfoil (Figure 1) are also very interesting in subsonic flow. 
During a very detailed study of dynamic qtall, the Ames USAAMRDL built 
a NACA 0012 airfoil with- a-’ 1.20 meter chcjrd for its 7 x 10ft wind tunnel 
(References 6 and 7). The tests covered ^Reynolds numbers ranging from 

g 

1 to 3.5 x 10 for the NACA 0012 airfoil and for other airfoils obtained 
by changing the shape of the leading edge of that same airfoil. 

A movie to be presented ' during this -symposium shows: 1 '.-':* 

- a case of leading edge dynamic stall of the type called "bursting of 
the leading edge laminar bubble when the radius of the leading edge of 
the NACA 0012 airfoil is reduced to a 0.4% instead of the standard 1,58, 
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- a case of trailing edge stall of the type called "gradual flow 

reversal of trailing edge separations toward the leading edge" 

when the NACA 0012 airfoil is equipped with a cambered extension 
designed by ONERA, I 

- a case of complex stall for the reference NACA 0012 airfoil as 
it involves sudden turbulent leading< edge stall preceded by a 
non-negligible reversed flow of trailing edge separations. 

We will cover in detail these different types of dynamic stall in 
the next paragraph but it is to be noted that in each of these cases, 
we have, in the critical phase of generalized separation, the presence 
of vortices that are more or less organized and complex as shown in 
the photos of Figure 6. 



Figure 6 -Evolution of separations of a vortex nature with an oscillating 
airfdil . 

a) NACA 0012 with a reduced leading edge radius 

b) NACA 0012 with a cambered extension 

Re B 2.5 x 10 ;a=15°+10°sin wt; k=0.15; Ames USAAMRDL tests, 
c 
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2.2 - Unsteady- Separations of the Turbulent Boundary Layer. 

Heated films are good detectors of heat flux. The convedted flux 
depends very -much -on the laminar , ..turbulent or separated condition 

t 

of the boundary layer at a perpendicular angle to the heated .film 
(reference 8) . 

Thus, as shown in figure 7> during [S 10 wind tunnel tests at the ToulouseJ 
Aeronautical Test Center, we showed a gradual flow reversal of the 
turbulent separation originating from -the NACA 0012 trailing edge Gaith 
a cambered extension to the airfoil) oscillating in pitch on both sides 
of its steady stall angle of attack .(Reference 9). This phase precedes 
the generalized separation phase with vortices as described in th'e 
preceding , paragraph . 

The same type of information can be obtained with heated wires placed 
very close to the top of the airfoil. This technique has been used by 
the Ames USAAMKDL (Reference 10) during the test program- with the 
NACA 0012 airfoil and its derivatives, a test program already mentioned 
above. Figure 8, extracted from reference 2, shows the development of 
separations from various airfoils tested against time or the geometric 
incidence angle of the airfoils. To obtain this figure } we have made 
the point where separation takes place the same as the point where 
reverse f low.sappears for the first time, a procedure that is not very 
rigorous. 

However, it is clear that, under identical' conditions of oscillations 
in pitch: 

- the NACA 0012 airfoil with a reduced leading- edge 'sees separations 
take shape and' propagate maiply from the leading edge to the trailing 
edge, 

- the NACA 0012 airfoil with a cambered extension sees separations take 
shape and propagate from the trailing edge to the leading edge. 





Figure 
with a 


Figure 





7 - Evolution of turbulent separations with a NACA 0012 airfoil 

g 

cambered extension. M=0.2; k=0.114: Re =1.9 x 10 . 

r 


~ TIME OH ' 7 

FLOW 

REVERSAL, wt 

' 7r/2 £ \ CLASSICAL TRAILING . - 
\ / EDGE STALL 

X (ONERA CAMBER) 
ooig- 


INCIDENCE, 

a, deg 


ir/4 




— rr/8 \- 


0 012 TR1R . 

'ABRUPT TURBULENT 
SEPARATION 



| ^ 

LEADING EDGE \ 
BUBBLE BURSTING 
(SHARP LEADING EDGE) 

1 i i . 


^25 • Opv 

24 . 

' , / 

22 V 
20 

18 ' ’’ 
15 « a Q 
12 


-Ca-j 

1.0 


8 - 


o - .25 .50 .75 

CHORDWISE LOCATION. X/C 

Evolution of turbulent separations with ah “oscillating airfoil. 

a=15°+10°sinu>t; k=0.15; Re=2.5xl0 . 


-10 



- the NACA 0012 airfoil first sees a trailing edge separation move back 
toward the leading edge then sees a sudden separation take place from 
the leading edge before trailing edge separations can reach the leading 
edge whether the transition has been triggered or not. 

These separations of the turbulent boundary layer depend, as in stationary 
flow, on many parameters: 

- those that condition the value of the local pressure gradient: 

i 

- the airfoil: itself as illustrated in the preceding Figure 8, 

- the average incidence angle and the amplitude of oscillations in 
pitch, 

- the reduced oscillation frequency since the time derivatives of 

V 

the airfoil incidence angle are directly proportional to it. 

f 

Figure 9, derived from Figure 26 of reference 10, illustrates the effect 
of this parameter. The higher the reduced oscillation frequency, the 
later separations appear. This-delays the generalized •■separation over 
the airfoil by the same amount. 


“ those associated with viscosity: 


- the Reynolds number 


- the location of the boundary layer transition, whether ^created 


naturally or artificially. 
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Figure 9 - Effect of the reduced frequency on the evolution of separations. 
NACA 0012^ a=15°+10°sinajt 
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Figure 10, derived from Figure Jt xn reierente xu illustrates the 
effect of the Reynolds number. ' The higher that number, the laterr 
the separations appear and develop. 

The effect of upstream conditions on separation will be discussed in 
Paragraph 4,3. 


cot <&*, inciDEnct ’ 

War a iw = 



Figure 10 - Effect of the' Reynolds number on the evolution of separations. 

NACA 0012. a=15°+10°sinajt ; k=0.15. 

2C - .abscissa of the separation point. 

t 

3 - TOTAL FORCES' AMD LOCAL PRESSURES ON AIRFOILS IN A DYNAMIC STALL. 

I 

3.1. - Typical Example of an Airfoil with- Harmonic Oscillations in. Pitch. 

Figure 11, extracted from reference 2, picks up the lift and pitching 
moment cycles of NACA 0012 airfoil-oscillating in pitch already depicted 
in Figure 1 but it links the main consequences s * ,of dynamic stall 
described/ in detail in the preceding paragraph to different times in 

I i 

the incidence cycle. 

We will particularly note that the lift of the airfoil continues to 
increase with the incidence angle, way beyond the maximum value obtained 
in stationary flow and this in spite of the presence of reversed flow 
inside an 'increasingly turbulent and thick boundary layer over a large 
section of the chord of the airfoil (a). A' vortex is generated near 
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(o) BOUNDARY LAYER FLOW 
REVERSES, LARGE 
lEDDIESI -APPEAR 

(d) MAXIMUM NEGATIVE 
MOMENT 

(b) VORTEX FORMS NEAR 
•LEADING EDGE . 

(e) FULLY SEPARATED 
FLOW 

(c) MAXIMUM NORMAL 
FORCE 

Hi . ' - 1 

(f) BOUNDARY LAYER 
REATTACHES FRONT 
'TO REAR 


Figure 11 - Main events associated with the NACA 0012 airfoil in a 

g 

dynamic stall .configuration. Re = 2.5 x 10 . 

the leading edge (b) and the lift (.-continues to increase with increased 
speed with a greater.incidence angle. The departure of the vortex toward - 
the trailing edge causes a moment stall of the airfoil but the lift 
stall . Cc intakes place only a little later when the vortex has gone, 
in this particular case, a little beyond the half-chord of the airfoil. 
Negative pitching moments (Cm < 0) can take take much higher values than 
during steady flow stalls and this is due primarily to loads induced by 
the vortex during its motion. The maximum negative moment corresponds 
to the passage of the vortex over the trailing edge (e). The incidence 
angle continues to decrease, the boundary layer finally reattaches 
itself to the airfoil, starting with the leading edge (f) but we find 
a non-separated airfoil only at incidence angles much lower than the 
incidence angles for stall in a steady flow. 
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Absolute pressure measurements made during tests "conducted by the 

. * \ 

Ames USAAMRDL are very useful 'in following the evolution of phenomena 
associated with the presence, development and motion of vortices that 
are characteristic of dynamic stall. 

Thus, Figures 12 and 13 extracted from aj document to be published 1 
(reference 11) show that: 

- pitch stall (b) is caused by the sudden reaction of high pressure 
drops at the leading edge (a) (-Figure 12) but is slightly preceded by 
that drop. Ther pressure dfop takes pla<j:e earlier but less suddenly 
if we trigger the transition artificially. Lift stall (c) takes place 
only a little later. 
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Figure 12 - Evolution of pressures in the leading edge area. Case of 
the NACA 0012 airfoil in a free transition, 
a = 15° +10°sinmt; k=0.15. 
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“ the vortex moves along the topside, of the airfoil generating local 
pressure drops as it moves (Figure 13); the maximum pressure drops make 
it possible to locate the vortex; we find that the displacement velocity 
of this vortex is generally less than half the flow velocity at an 
infinite upstream distance. A detailed study of this phenomenon is also 
made in reference 12 by F.O. Carta, using tests made by Sikorsky in 
1972 using a NACA 0012 airfoil. 
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Figure 13 - Effect of the pasfeing-'vortex on pressure distributions. 

/ 

Case of the NACA 0012 airfoil with a fixed transition, 
k = 0. 15; a - 15° + ,'l0° sinait. 


Knowing the local absolute pressures makes it also possible to follow 
the evolution of pressure drags. Figure 14 shows the pressure !drag 
cycle along the axes associated with the model, ' 


The absolute pressure measurement technique has also been used by other 
researchers such as A.G. Parker (reference 13) and previously by R.I. 
Windsor (reference 14); these works are very useful as they also provide 
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Figure 14 - Effect ^pf the reduced frequency. 

airfoil with a free transition. 
a - 15°+ 10° sin ai t. 


Case of the NACA 0012 


a basis for comparison for developing computation' methods for unsteady 
flows in 'the presence or absence of dynamic stall. 

3.2. - Characteristic aerodynamic value' for -airfoils with harmonic 
oscillations in pitch. 

The typical previous example reveals a certain number of aerodynamic 
values for which it is interesting to know the evolution as a function 
of oscillation parameters." The - following are mainly involved: 
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“ the maximum lift reached by the airfoil' during an oscillation since 
it quantitatively reveals the beneficial consequence of the delay in the 
stall as noted experimentally, 

- the maximum negative moment obtained during th'e critical phase of 
dynamic stall since it can quantitatively reveal the penalty level 
of dynamic stall. 

Another significant point of dynamic stall configurations involves 
energy transfers between the "air foil and the airflow which, therefore, 
condition aerodynamic stability. 


The inside area of the pitching moment cycles as a function Wf incidence 


angles is proportional to the work of the aerodynamic forces (f C TO da) 
If the cycle is accomplished in a counter-clockwise direction, -the- 
aerodynamic damping is positive and the motion is stable. On the 
^other hand, if the cycle is accomplished in a clockwise direction, the 
work by the aerodynamic forces is positive, the aerodynamic damping is 
negative and this can lead to a tension instability of the blade. 


The very complete tests conducted in 1968 and 1969 by Boeing-Vertol 
with 4 airfoils (references 16 and 17) equipped with differential 
pressure sensors and moving in pitch and up and down with a varying 
reduced frequency showed that the amplitude and the average incidence 
of the oscillations and the Mach- number make it possible to determine 
the effect of these parameters on the unsteady aerodynamic performances 
of the airfoils. 


There is no question here of providing all of the, typical results . 
obtained .but simply to bring tic' attention to the primary effects, 
particularly the followings- 

- th&- average incidence of the oscillation and the associated amplitude 
as the effects are very different depending on whether the range of 
angles examined is located ahead of, on both sides or behind the steady 
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flow stall incidence angle of the airfoil. Figure 15 shows an example 
of the cycles obtained. We will note that, it is in an oscillation on 
both sides of the steady flow stall incidence angle that the effects 
are the most noticeable and that they lead to large maximum lifts and 
to the possiblity of instability in pitch. When the airfoil oscillates 
continuously in separated -flow, with high average incidence angles, 
motion 'becomes stable again; however, generally speaking, this takes 
place with high values of negative, pitching moments. 



original page ft 
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Figure 15 - Cycles of the forces and moments obtained with a Vertol 

g 

23010 airfoil with harmonic oscillations in pitch (Re^ = 4.8 x 10 ). 

- the Mach number, since as soon as this parameter- increases a little 
and we have high incidence angles, problems associated with compressibility 
appear. Figure 16, extracted from reference 2,- shows the cycles obtained 
with different Mach numbers and' with average incidence angles that are 
slightly greater than the corresponding stall incidence angles 'for steady 
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Figure 16 - Effect of -the Mach number on the unsteady flow characteristics 
of the NACA 0012 airfoil. 


flow. If the general shape of the cycles^ remains the same, we will note 

that the maximum lift decreases with the Mach number and that', in general - , 

i 

the maximum negative moments are less significant at a - ’ high Mach number. 

We will note particularly the tendency of pitch instabilities at a high 
Mach number to disappear. This can be caused by a weaker generation of 
vortices when the shock wave is^ located near the leading edge of the 

j 

airfoil. ' 


These tests have been analyzed by many authors (references 2, 18, 19 and 
20 being given as examples) . 

The effect of the shape of the airfoils tested is not. negligible;; this 
is illustrated in Figure 17 which regroups Figures 7, 8 and 9 of reference 
19. 


We will note that: 

- within the range of the thicknesses tested, the lower the airfoil thick- 
ness, the lower its maximum lift capabilities in unsteady" 1 . flow (as well' as 
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STUDY OF BLADE STABILITY j 
IN PITCH 



Figure 17 — Effect of the thickness and camber of the airfoil. . Boeing- 
Vertol tests with airfoils oscillating in pitch, a - _+ 5° 

i 

in steady flow). Moreover, the average incidence angle at' which the 
airfoil becomes unstable decreases with a decrease in airfoil thickness. 

The instability pocket of thin airfoils extends closer to the high 
incident Mach numbers than is the case for thick airfoils. 

- the effect of the camber is generally favorable and, as is the case 
in steady flow, this effect provides greater maximum lifts than are 
available with symmetrical airfoils with an equivalent thickness. Pitch 
instabilities then appear at higher average incidence angles for a fixed 
Mach number, qr at a higher Mach number with a fixed average incidence angle. 
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We have been able to observe an ’’instability pocket” in the average 
incidence/Mach number plot for the test. This involved oscillations 
with average, amplitudes (+5°). Phenomena associated with both compressibility 
and stall being highly non-linear, we must expect these pockets to depend 
on the amplitude and pitch instabilities to exist at Mach numbers 
greater than Mach 0.6, for example, 'with low-amplitude oscillations. 

Figure 18 gives an axample of this and is extracted from the stall flutter 
test results obtained from tests conducted by the Direction des Structures 
(Structures Management Office) of ONERA (reference 21). 





Figure 18 - Example of a result obtpinejd in a stall flutter study. 

Oscillation axis located at -0.375c. j 
NACA 63A015 oscillating, in pitch - a^^lO 0 . 


Outside of the tests already mentioned, we can obtain information on 
unsteady flow characteristics of other airfoils' thanks to tests conducted 
by the RAE on -the NPL 9615 airfoil for example (reference 22) or by the 
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■ Toulouse CEAT (Toulouse Aeronautical Test Center) in France using the 
SA 13109 airfoil for example (reference 23). New test benches designed 
more specifically for high Mach number tests are being assembled or are 
just entering into operation. These include, as examples, facilities 
planned at Ames and Langley Field by NASA, or at the ARA of Bedford 
under the RAE. 

3.3. - Other motions. 

If the analysis of dynamic stall phenomena was performed on changes 
of model incidence angles, it is because we have a lot of information 
available on the motion associated with these Wangle changes. 

We note also: 

- the up-and-down motions of the airfoil] obtained by subjecting it to 
cyclical variations in translation perpendicular to the upstream flow 
velocity at infinity IT . A few typical characteristics of the tests 
performed by Boeing-Vertol (references 16 and 17) are provided in 
references 18 and 19. However, we will note that this motion is 
equivalent -as a first approximation- to a variation in the airfoil 
incidence angle with a hto/U^ modulus, h being the amplitude of the 
up-and-down motion and w being its pulse.. Differences are caused by 
the fact that, over the airfoil, the distribution rule for the 
velocities induced, by the motion is linear along the 1 pitch oscillation 
motion whereas it is constant along the up-and-down motion. Figure '19 
gives an example of the cycles obtained and confirms the possibility 
of obtaining unsteady lifts thatf are still greater than steady flow 

I i 

stall lifts but at the costs- of large negative moments and of motion 

instabilities. Figure 20 depicts the evolution of the primary values 1 

that are characteristic of dynamic flow configurations as a funtion of 

/ 

t 

the airfoil setting angle and of the upstream Mach number at infinity, 

- the back-and-forth motion of the airfoil. The motion is .then a cyclical 
translation of the airfoil in a-'direction parallel to the upstream air 
flow velocity at infinity U^. In fact, we obtain a cyclical modulation 
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Figure 19 - Example of unsteady cycles obtained by rtioving the Vertol 
23010-1.58 airfoil up and down. 



23010 airfoil with up-and-down oscillations. 
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of the airfoil attack velocity reproducing, for example, the conditions 
of sinusoidal attack velocity for a helicopter blade section during a 
translation flight. It is the Marseille IMF that is working in this 
area in France (references 24 and 25). There again, there is a possibility 
of obtaining, during these motions; lifts that are greater than those 
obtained with steady flow, depending on the incidence angles to which 
the airfoil has been set (greater or equal to the stall'indidence angle 
under steady flow conditions), the reduced frequencyyof the motion 
k^ffc/U^ and of its reduced, amplitude A- Aw/U^A being the actual amplitude 
of the back-and-forth motion). Figure 21, derived from Figure 1 in Reference 
26 illustrates this point but the paper presented by the Marseille IMF 
.engineers during this meeting (reference 27) will make it possible to 
cover in greater detail the flows associated with this motion and the 
aerodynamic performances thus obtained. This attack velocity variation 
can also take place by varying the flow velocity inside the tunnel with 
the airfoil remaining fixed. This' type of installation has been adopted 
in the United States at the University of California at Monterey, at 
the University of Texas and at the University of Illinois. Tests performed 
at the latter university using a NACA 0012 airfoil will be presented 
during this meeting by L.S.Saxena, A.A. Feher and M.V.Morkevin (paper 



Figure 21 - Back-and-forth harmonic oscillations with a NACA 0012 airfoil. 

1 < Re < 2xl0' > . 
c 
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4 - PREDICTION MFTHOPS 

To predict the complex phenomenon that dynamic stall represents, we have 
different methods available. Some bring only partial information; 
others pertain to potential flows with or without displays'‘of unsteady 
boundary layers and separations. The others, based on the solution 
to the Navier-Stokes equations are promising but are still in limited 
use. We also have available practical methods of an empirical nature 
when we are only interested in overall unsteady flow forces. 

This paragraph will attempt to provide sjome information on the various 
available methods illustrated through several characteristic results. 

4.1. - Computations of unsteady potential flows in the absence of separations. 

When airfoils are not subjected to large viscous effects or separations, 

i 

the solution of unsteady equations for perfect fluid mechanics makes it 
possible to describe in an adequate manner the pressures, forces and 
moments applied to the airfoils. 

Basic theories apply to airfoils with zero thickness and harmonic 
motion. Theodorsen 1 s’ '‘■theory (Figure 15), for example, reveals naturally, 
in its formulation, the notion of motion reduced frequency (k-wc/210 and 
its effect on the modulus and phase of the pressures, forces and moments 
relative to the imposed motion. However,- only the terms corresponding 
to the fundamental harmonic of gj are computed and they are proportional 
to the amplitude of the motion. 

A simple extension qf this method to the case of thick airfoils has been 
proposed by W.J. McCroslcey (reference 2§) who provides velocity distributions 
on top and below oscillating airfoils. Associated with the solution to the 
Bernoulli equation, the method shows the existence of w terms which then 
depend on the shape of the airfoil (camber and thickness effect) , of " the 
average incidence angle, of- the amplitude of the reduced frequency 
and also of the amplitude of the 2w terms which depend on the amplitude 
and reduced frequency of the motion. This method has provided good 
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matches with test "results obtained by R.I. Windsor (reference 14) and 
can be used to perform computations for unsteady boundary layers, for 
example. 


The advantage of thes'e methods is to provide a mathematical formula 
that is relatively simple once the Theodorsen functions are known, 
these fimefiens depending ©n the only reduced frequency. 


Geising (reference 29) has perfected a much more sophisticated computation 
process by applying a method making it possible to handle any type 

j 

of airfoil in any type of unsteady motion. Computations can only be 
performed on a computer, starting with an airfoil at the rest position. 

At every time period, we compute the intensity of the sources positioned 
on the airfoil and of the associated vortices located inside the airfoil, 
and the position of the free vortices released from the trailing edge 
to form the unsteady wake, the development of which is also computed 
in space. Taking the Kutta-Joukowski condition, into consideration, 
this computation makes it possible to predict the evolution of pressures, 
forces and moments acting on the airfoil. If a few computation procedures 
for unsteady flow ‘in' the _ absence of any separation from the airfoil have 
been recalled here, it is because-:. 

- for any dynamic stall configuration, there are times when there are 
no separations from the airfoil, 

- the analysis of the results provided by these methods can partially 
explain stall delay. 


The Theodorsen theory shows through an example that when the incidence 

angles of an airfoil, oscillating in pitch increase, the .pressure gradients 

along the- chord are less than the pressure gradients obtained in a steady 

flow. , : Can this observation made 'by F.O. Carta -(reference 30) be sufficient . 

to explain the stall delay? Figure 22, .extracted from reference 1, shows 

the result obtained by assuming that the unsteadysstall (a ) takes place 

• JJb 

only when the unsteady pressure gradient near the airfoil leading; edge 
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(at x/c=0.1) reaches the theoretical value [for steady flow_ separation. 
Computations properly predict, as is actually experienced, that the 
unsteady stall incidence angle increases with the reduced frequency 
of the motion but computations are far from being able to predict the 
absolute value of the increase in the stall incidence angle in 
dynamic stall configurations. 



origjnae E*3B jj» 
pe eook oauJ® 


Figure 22 - Increase of the unsteady stall incidence angle a against 

* Ub 

the steady stall incidence angle a . 

bb 


4.2. - Computations of dynamic stall configurations using potential flows. 


We have "seen that, experimentally, one of the major phases of stall was 
the formation and displacement of a vortex on the top side of the 
airfoils. 

- Professor Ham (reference 31) has simulated this phase by generating 
near the -leading edge a series of vortices and by letting them escape 
toward the trailing edge at the local flow velocity. This- computation 
has been made for the case of an infinitely thin airfoil. 
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Figure 23, extracted from reference 31, shows the trajectories followed 
by the vortices in the case of a flat plate suddenly subjected to an 
incidence angle of 30°. We can readily see a vortex mass being 
generated and moving toward the trailing edge. 
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Figure 23 - Computed trajectories for the vortices generated when the 
incidence angle is suddenly set to 30°. 

- i 

Figure 24, also extracted from reference 31, shows a satisfactory 
agreement with actual experience with, respect to local pressures on 
an airfoil oscillating in pitch between 10 and 30° at a reduced frequency 
of 0.35 (tests performed at MIT). 


However, it must be noted that we must assume the time at which leading 
edge vortices are emitted. Moreover, computations cannot .take into 
consideration trailing edge separations. 
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Figure 24 - Computed and experimental evolution of the pressures on 
an airfoil in a dynamic stall configuration. 


- In' France, the Societe Bertin (references 32 and 33) has used this 
scheme hy applying it to the case of [an jairfoil. But it is the separation 
of the boundary layer ’ taking place toward the trailing edge that 
generates a vortex layer. This layer is generated from this point and 
is fed with a rotational flow "by. the boundary • layer . According to 
Sarpkaya (reference 34) , the rotational flow, rate is given hy the 
formula: • Ay/At = -1/2 (U^-U^) . 
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where U is the flow velocity at a right angle to the point of laminar 
separation and is the maximum velocity of the recirculation flow 
induced by the vortex winding action. 

At every point in time, we therefore emit discrete vortices for which 
we compute the trajectory. We determine the velocities induced at the 
airfoil wall by all of the singularities corresponding to the J.P.' Geising 
method (reference 29) and by the free vortices generated at the leading 
•'edge . 


Figure 25, extracted from reference 33, shows a comparison between computed 
vortex lines and smoke visualizations performed by the Ames USAAMRDL 
(reference 6).' Comparisons .of lift cycles obtained are rather encouraging, 
as shown in Figure 26 extracted from reference 2. 




Figure 25 - Evolution of uns_teadg separation. Comparison of computations 
and visualization, a = 15°; a- 6°; k= oi24i 
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However, there remain" some problems in determining either rigorously 
or empirically (under study) the conditions for the generation and 
diffusion of these vortices. 

Potential flow computations still have their own limitations and only 
unsteady boundary layer computations canj bring information elements to 
make it possible to understand and possibly predict the delays: in . 
the appearance of separations. 



Figure 26 - Comparison of Bertin computations and USAAMRDL (Ames) '".tests . 
NACA 0012 airfoil. 


4.3. - Computations of unsteady boundary. layers . 

One of the sessions in this unsteady aerodynamics meeting is devoted to 
the problems of the unsteady boundary layer, whether separated or not; 
it is preceded by an overall presentation by D.P. Telionis; therefore, 
it is not useful for us to develop this subject' here. 
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Here, we will only recall a few fundamental problems associated with 
dynamic stall configurations. 


Figure 27, extracted from reference 1,. shows that in the case of the 

oscillating airfoil, the effect of the unsteady terms 3u in the turbulent 

9t 

boundary layer equations is not negligible in the area where the drag 

coefficient becomes zero. In the computations that were performed, 

we took the velocity distribution laws fjor potential flow as determined 

by W.J. McCroskey (reference 28). The hysteresis effect which is observed 

' 3 u 

to be more significant when we consider the ~ terms comes primarily. 

d t 

as the result of the fact that - near the 'trailing edge, these time derivatives 
are of the same order of magnitude as the gradients along the chord 


3u 
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Figure 27 - Evolution of ,the point of zero drag -coefficient with an 
oscillating NA.CA 0012 airfoil. Re-3xl0^; k=0.24; a-15°+ 6° 
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A detailed analysis of the reversal toward the leading edge of the 

point where 1 a reverse flow appears (equivalent to the point where the 

drag coefficient is zero) was performed by R.M. Scruggs, J.F. Nash and 

R.E. Singleton (reference 35). Figure 28, extracted from this last 

reference, shows the computation results obtained with an airfoil 

subjected to an increasing incidence angle variation performed at. 

a constant velocity and between 0° and 25°. If we consider the computations 

performed with the unsteady potential flow and the unsteady boundary 

limit, the point where the reverse flow appears starts moving toward the 

■ -- . ' 
leading edge later but, with an increasing incidence angle, it moves much 

faster than in steady flow. The delay in this motion is caused, according 

to the authors, in part by the unsteady 'potential flow and in part by 

unsteady boundary layers as shown in the different curves of Figure 28. 







Figure 


nation ox the point 


(Increasing incidence angle and ca/u^Q .04) 
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It is significant to note (Figure 29, extracted from reference 35) 
that if .\ze consider the incidence angle at which the point where 
reverse flow takes place appears at mid-chord, this angle will increase 
with the speed of airfoil .angular variation in a manner very similar to 
lift with respect to the steady stall incidence angle.' This increase 
is at first very steep for low values of ca/211^ and then there appears 

a type of saturation of the unsteady effects at high values of cct/21 J^. 

This characteristic phenomena' reportedly originates also in part from 
the unsteady potential flow and in part from the unsteady boundary layers. 



Figure 29.- Incidence angle at which reverse flow reaches the airfoil 
mid-chord as a function of the angular velocity of the airfoil. 


- 34 - 



However, we must state that these computations can really be meaningful 
only if the point where reverse flow appears is not followed downstream 
by too large a boundary layer separation which would then change velocity 
distributions' over the entire airfoil (this’ would appear to be the case 
with tests performed at Ames using the NACA 0012 airfoil with or without 
its cambered extension), 

A major point associated with the unsteady nature of boundary layers is 
the distinction between the' point where reverse flow appears and the 
point of separation. Very many studies have been conducted in this 
subject area and the reader can refer to certain papers given as examples: 

- papers presented at the IUTAM Symposium on unsteady boundary layers 
organized by Laval University, Quebec, Canada, in May 1971 (reference 36), 

- papers presented at the symposium organized by the University of Arizona, , 
Tucson, in March 1975 (reference 37) , on unsteady aerodynamics. 

In fact there does not appear to be any singularity, at the point of 
reversed flow as in .stationary flow. For exampley the separation criterion 
is involved only upon the appearance of. the divergence of the 'displacement 
thickness. This takes place upstream from the point where reversed flow 
appears as shown in the sketch of figure; 30 proposed by W.R. Sears and 
D.P. Telionis (reference 38). Therefore;, we can continue 'to describe 
the boundary layer upstream of .the point of reversedflow (cf=0). V.C. 

Patel and J.F. Nash were able to compute (reference 39) the change in 
time of a turbulent separation bubble created by a deceleration followed 
by an acceleration of the external flow which is variable in time. Figure 31, 
extracted from reference 2, illustrates the results obtained which show 
the birth of the separation bubble with its reversedflow, its extension 
toward the leading edge and the final appearance of separation caused 

V 

by the sudden increase of the boundary layer thickness. These phenomena 
of unsteady turbulent separation bubbles are also studied in an intensive 
manner by L.W. Carr of. the Ames USAAMRDL. 
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Figure 30 - Sears and Telionis model for the separation of a boundary 
layer moving back toward the leading edge. 



Figure 31 - Computation of the variation in the displacement thickness 
of an unsteady turbulent boundary layer with reversed flows. 
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There are many methods for computing unsteady turbulent boundary layers 
but we must mention that the selected turbulence models are, for the 
time being, always based on criteria that are valid ' for steady flow. 

It would therefore appear fundamental to this -field to check that these 
models are suitable to unsteady aerodynamics. Experiments based on 
research of unsteady boundary layer characteristics will therefore be 
very useful and we will mention, as an example, the work performed by 
the Toulouse CERT-ONF.RA with flat plates subjected to pulsed flows. This 
work will be the subject of a presentation during session III of this 
meeting (reference 40). 

In conclusion, boundary layer and potential flow computations provide 
interesting partial data and make it possible to understand certain 
phases of unsteady stall. However, they remain incomplete. 

4.4. - Navier-Stokes Equations Computations. 

In a problem as complex as the interaction between potential and viscous 
flows, the solution of the Navier-Stokes equations should make it possible 
to approach the dynamic stall problem in its entirety. 

In solving the vortex propagation equation and Poisson's equation'‘for the 

stream function through a scheme using finite differences and a relaxation 

procedure, U.B. Mehta and Z. Laven had already used the stall of an airfoil 

after it had suddenly been set to an angle of attack (reference 41). A 

similar approach' was used for the case of an oscillating airfoil and 

U.B. Mehta will present, during this symposium, the results that he has 

obtained (reference 42). Visualizations obtained at the ONERA aerodynamic 

tunnel under conditions corresponding to those used in the computations 

make it possible to dete.rmine the accuracy of the preliminary results 

obtained. Figure 32 shows a very good agreement between the computed 

’ 1 

streamlines and the paths of the air bubbles used to make visualizations 
possible. The vortex phenomena are properly predicted and we can therefore 
hope that computations will be able to describe in a satisfactory manner 
the physical process involving the birth, convection and eventual dissipation 
of these vortices' which are characteristic of unsteady flow. 
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It is hoped that computations can be made • : with more realistic Reynolds 
numbers by using turbulence models that are suitable to cases involving 
strong viscous interactions. 





Figure 32 - NACA 0012 oscillating in pitch. 

a =10°+10°sin wt; k=0.25; RE = 10,000 

0 

During this meeting, other papers also cover the solution of the Navier- 
Stokes equations. The paper by J.C. Wujand N.E. Sankar also involve 
an oscillating wing (reference 43) and the paper by R.B. Kinney covers 
the case of a fixed wing subjected to uhsteady flows (reference 44). 

The methods used to solve the Navier-Stokes equations may be different 
and it will be interesting to learn if some methods can be applied to 
higher Reynolds numbers without requiring prohibitive computation times. 

Another application of the Navier-Stokes equations was made by S.J. 
Shamroth and J.P. Kreskowsky (reference 45) by using them only near the 
leading edge. Computations provide a rather detailed description of 
the leading edge laminar separations and could therefore be used for 
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leading edge dynamic stalls through, for I example, the bursting of 
separation bubbles. 


4.5. - Problems caused by the coupling between viscous and inviscid 
flow. 

We can also use the so-called coupling methods which make it possible 
to compute pressure distributions by performing perfect fluid computations 
taking into consideration the development of viscous layers through the 
use of conditions at the limits. 

Figure 33 illustrates a trailing edge turbulent separation configuration 
with reversal toward the airfoil leading edge. In one coupling method, 
it will be necessary to take into account the laminar boundary layers, 
the transition, the turbulent boundary layers, the separation areas and 
the wake. 

For example, to take the boundary layer into account different methods 
already applied to stationary flow are possible: 

- we can redo potential flow computations after enlarging the airfoil 

with the computed displacement thickness using a boundary layer program 

and then iterating until there is convergence; 

'we could also, distribute along the airfoil a series of sources and "wells" 

corresponding to the flow rate variation affected by the boundary layer; 

the intensity of these singularities (q=d/ds(p U 6 ) is computed from 
■ e e 

computation results for viscous and potential flows in the preceding 
iteration. 

The transition takes place in many cases by means of a' laminar separation 
bubble. Unsteady conditions to which the airfoil may be exposed can 
retard its bursting. Work performed by J.D. Lang (reference 46) is, 
for example, interesting in this regard. The purpose of this separation 
bubble can certainly become important in the case of airfoils which, 
in steady flow, stall through this, bubble bursting. We must also mention 
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Figure ' 33 '~~ Exa.mole trf -viscous areas over an airfoil in a dynamic stall 
configuration. 

the studies in this area performed at MIT by N.D. Ham and his team 
(references 72, 73, and 76). ' 


The problem of the separation area and of the wake remains much more 
complex and deserves in-depth study. We can mention the studies by 
P. Crimi (references 47 and 48) who, in his work for predicting forces 
and moments on airfoil in a dynamic stall configuration, makes the 
hypothesis, for example, of a constant pressure in the separation bubble. 


The separation point on the airfoil is assumed to be known but is 
redefined at every consecutive iteration between computations of the 
viscous and potential flows until there is convergence at that point. 


by assigning a law .for the'-displacemen’t of the separation point on an 
oscillating airfoil, F. Sisto and P.V.K. Perumal (reference 49) have 
computed the evolution of the forces and moments. For potential flow 
computations, they used a potential acceleration method. Computations 
made for a flat plate oscillating in pitch and moving in a pumping 
action show that we can thus obtain force and moment cycles that can 
lead to flutter. In fact, the selection of singularities to depict 
separation areas and the viscous wake coming from the merging of the 
boundary, layers on the upper and lower surfaces seems to demand further 
in-depth analyses. 
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Nevertheless,, the theoretical aspect of this problem has been recently 
examined by W.R.' Sears (reference 50) who covers in particular the 
delicate problem of the Kutta-Joukowski condition in unsteady flow. 

It would perhaps be interesting to develop methods similar to those 
developed for predicting the maximum lift of airfoils in steady flow as 
B. Maskew and F.A. Dvorak (reference 51). have just done. 

The problem of coupling between viscous and inviscid flows will certainly 
require- much effort but the goal to reach is most attractive' since it 
should provide a much more accurate computation of dynamic stall configurations 
and of the aerodynamic qualities of the airfoils under these conditions. 

4.6. - Empirical Methods for Predicting Overall Forces. 

The preceding paragraphs show the different fundamental problems that have 
to be resolved to describe in an 'accurate manner the various phases of 
dynamic stall. For years, with a goal of using them for applications, 
many researchers have devoted themselves- to finding simpler methods for 
representing unsteady forces with airfoils . encountering dynamic stall 
regimes during unsteady motions . 

The problem consists primarily of finding a very general formula from 
tests performed with simple configurations (oscillating airfoil, for 
example) while varying many parameters: airfoil, Mach Number, incidence 
angle and amplitude, etc., 

UARL Method (a. A, and B) . 

This method has been perfected by F.O. Carta and rests on the fact that 
instantaneous lift and pitching moment must be functions -as in non-separated 
unsteady flow- of the incidence angle a, of the angular speed d and of the 

i»i» 

angular acceleration a. 

In fact. Carta uses, for greater generality, the following reduced 

2 2 • • 

parameters: o=a/a ; A=dc/2U' ; B=6£c /4U where a ^steady stall .incidence ang. 
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We have tables available (references 52 and 53) which provide for 
any incidence angle average or C^ values for corresponding A and B 
values . 


This method has made it possible to resolve many problems associated with 
dynamic stall configurations or with configurations where aerodynamic 
instabilities appear (flutter) on a compressor or helicopter rotor 
(references 54 through 56). 


More recently, R'.L. .Bielawa from the United Technologies Research Center 


in Hartford (reference 57) has developed a method that makes it possible 
to break away from .tables using A and B values. Non-linear functions 


for which only the coefficients are tabulated synthesize 'the C^ 
curves = f(a, A, B). 


(or C ) 
m 


These methods can , for example, be applied to determine the angles of 
attack assumed by a model' totally' flexible in torsion when subjected 
to forced oscillations in pitch. Figure 34, extracted from reference 57, 
shows the good match that exists between experiments and computation results. 



Figure 34 - Comparison of the incidence angles assumed by a flexible 
model subjected to forced oscillations in pitch (<* 3^=5. 53). 
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Boeing-Vertol Method, 

This method was developed following a test program with airfoils oscillating 
in pitch and moving in a pumping motion performed by Boeing-Vertol 
(references 16 and 17).' It rests on the experimentally observed fact 
that unsteady effects are generally greater with an increase in the 
speed ip the ipcidenceapgle variation imposed on the airfoil. D.W. Gross 
and F.D. Harris (reference 58) have shown that the stall delay is 

r 

proportional to vcd/2U^. F.D. Harris and F.J. Tarzanin (references 59 
and GO) -then perfected a computation of aerodynamic forces by associating 
with the basic theory of Theodorsen a pseudo-dynamic or-'unsteady reference 
incidence angle which indicates the possibility of having stall. *• 
delay with an increasing incidence angle. 


R.E. Gormont (reference 61) has described in a very complete manner 
the Boeing-Vertol method which has the advantage, although based on a 
certain amount of empiricism, of a mathematical formulation which requires 
only that the airfoil steady poles be known. 


Figure 35, extracted from reference 60, shows that computations make it 
possible to predict with some accuracy the -cycles of forces and pitching 
moments for some airfoils in the absence or presence of total or partial 
stall. 



Figure 35 - Cycles of calculated and measured unsteady forces, 
(Vertol 23010 airfoil- M=o.4; k=0.‘062) 
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At.ONERA, this method has been modified to provide the capability to include 
in a program for computing local and overall performances of a helicopter 
rotor, the non-linear effects of unsteady stall configurations ■ (reference 62). 
Figure 36 shows that a more satisfactory agreement with changes in 
experimental local lifts has been obtained by using this method- (in comparison 
/with the purely linear original method) ‘on a rotor where a good portion 
of the retreating blade is in a dynamic stall condition. 


Dynamic Stall . Time-Delay Method,. 

This method is based on the simple idea that there could be two reduced 
time values and which determine the delay in moment stall and in 
lift stall. If t Q is the original time, at the time where the airfoil has 
an instantaneous incidence angle equal to the incidence angle where it stalls 
in steady flow, the times t^ and t 2 at which, under unsteady conditions, the 
airfoil stalls in moment and in lift, correspond to reduced times: 

U oo 


f N U CO \ 

V‘W x T 


WV * c 


where and are "universal" constants. 

To Sikorsky (reference 63) and have values of 2 and 6; to Westland 
(reference. . 64) , t^'and t have values of 2.44 and 5.41 after a detailed 
"examination of the Boeing-Vertol tests (references 16 and 17) for which 
these average values would be valid in more than 120 configurations. From 
this somewhat surprising result has come a simple method of computing 
unsteady aerodynamic forces in the case of dynamic stall. The four 
primary computation phases are illustrated in Figure 37, page 46. 


)As long as the airfoil does not reach the incidence angle corresponding 

1 

to unsteady moment stall, computations are made in potential flow without 

i 

separation. 

2) For intermediate times between moment stall and lift stall, we continue 
to compute lifts for non-separated flows butj the moment is computed by 
applying a linear law 1 as a function of the recoil time of the aerodynamic 
locus until we have a location equivalent to that obtained in totally 
separated steady flow. 
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3) when unsteady lift stall comes into play, we apply an exponentially- 
decreasing law of this lift as a function of time toward values for 
a quasi-steady pole of separated flow over an airfoil. 

!' 

A) finally, when the incidence angle oncje again becomes equal to or less 

i 

than the gtpady stall incidence angle, we return to the steady poles 
through a re-attachment process which is also expressed mathematically. 

With this method, x?e are also able to satisfactorily reproduce experimental 
lift and pitching moment cycles on airfoils in a stall configuration 
as shown in Figure 38 which has been extracted from works by T.S. Beddoes 
(reference 64). However, as in the case of the preceding methods, the 
agreement is not always good. 



Figure 37 - Primary phases of computations performed using the "Time- 
Belay" method. 
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Figure 38 - Effect of the reduced frequency of harmonic oscillations 
in pitch on the Vertol 13006 airfoil (M=0.2; (Xq= 12.5 0 ; a^=5°) 


Lockeed Method. 

L.E. Ericsson and J.P. Reding base their analyses of the results and 
their prediction method on a half-theoretical., half-empirical decomposition 
of the stall delay. 

Associated mathematical formulas are developed in reference 65. Their 

t 

application to the case of the tests performed at Ames by the USAAMRDL 
sometimes gives good results as shown in Figure 39 extracted from reference 
66. The phase for the appearance of moment stall is well predicted. The 
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phase for lift stall and the phase for the appearance of the maximum 
negative moment are associated with a vortex propagation speed of 
0.5511^ and to its passage over the 70% point on the chord or over the 
airfoil trailing edge. 
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Figure 39 - Prediction.<q.of events - associated with the dynamic stall of 
oscillating airfoils.. NACA 0012 airfoil oscillating in 
pitch. a=15°+10°sinwt. 


The methods developed by L.E. Ericksson and J.P. Reding have been the 
subject of many documents (references 67 through 70) bearing in particular 
on unsteady stall, flutter or dynamic effects associated with separation 
behind shock waves. Luring this meeting, these authors present a paper 
(reference 71) which states the problem of scaling up the results acquired 
in wind tunnels which must be transposed to full-scale aircraft 
within the scope of dynamic tests where the Reynolds number plays a very 
significant- role. 


MIT Method. 

It is based on the knowledge acquired at MIT by Professor Ham and his 

team in both the' experimental and theoretical areas (references 72 and 

73). W. Johnsson (references 74 and 75) has shown that the maximum lift 
* 

and the maximum pitching moment are linearly related to the angular speed 
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parameter ca/2U ro tut tnat there limiting values to their increase. The 
unsteady stall angle is estimated through semi-empirical formulas 
established by^-W. Johnsson and N.D. Ham |( reference 76). As an example, 
reference 75 provides the results obtained with this method and the 
ones developed at Boeing-Vertol and UARL in determining local loads in 
a helicopter rotor. ) 


W.J. McCroskey has provided, in reference 1, the results obtained with 
the various prediction methods mentioned above for determining maximum 
negative moments and lifts reached with 'the NACA 0012 airfoil tested at 
Ames with harmonic oscillations in pitch. Since the originial publication 
in March 1975, computations have been continued, particularly by Lockheed 
and MIT following improvements that they were able to make to their 
methods. Current results of the comparison between computations and 
experiments are provided in Figure 40. Overall, the results seem highly 
satisfactory but these are'very partial results. It certainly would be 
interesting to make a more complete comparison of the various methods 
by looking also at the phases during which the main effects 'on lift'" 
cyclesirand particularly the moment cycle take place. 



Figure 40 - Comparison of computations and experiments for normal forces 
and maximum negative moments reached with a NACA 0012 airfoil 
oscillating in pitch. 
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These various methods satisfy specific needs of manufacturers and 
researchers. The applications primarily involve helicopter rotors 
with certain flight configurations that impose stall conditions on . 
the retreating blade. Dynamic loads, stresses and deformations 
induced by the blade stalling are clearly much better predicted by 
using the preceding methods than was the case a few years ago. However, 
there still remain significant local or instantaneous differences which 
justify 'pursuing the search for new prediction methods that are increasingly 
sophisticated, or for unsteady tests with a three-dimensional effect, 
figure 41, extracted from reference 77, shows the test results recently 
obtained at Hartford with swept oscillating airfoils and the significant 
differences obtained with a 30° and a 0°. sweep angle, with respect to 
aerodynamic moment' cycles. It appears tjhat this wing sweep angle 
increases the tendency for the airfoil to be unstable in pitch. This 
illustrates the aerodynamic complexity that stall configurations must 
have in three-dimensional and unsteady flows. 
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5 - CONCLUSIONS 

Dynamic stall is a complex example of strong interaction between viscous 
and inviscid flows. Even in the simple case of oscillating airfoils, 
experiments involving very complete measuring and observation equipment are needed 
order to be able to properly interpret the complex phenomena associated 
with such eonfigurations , The physical process for forming and feeding 
vortices on top of the wing still remains poorly defined, hence difficult 
to show through computations. Means for predicting the evolution of unsteady 
forces and moments need to be perfected further. "Although more or less 
empirical formulas are at times satisfactory, it will be necessary to 
have programs for the coupling between viscous and inviscid flows which 
take into consideration the boundary layers, separations, the wake and 
the vortex phenomena that are characteristic of dynamic stall. 

Instabilities or flutters are one of the most important consequences of 
dynamic stall. Unsteady bench tests operating in transsonic wind tunnels 
will make it possible to have tests that are necessary to the proper 
understanding of the -effect of shocks and separations on these phenomena. 

The three-dimensional aspects of dynamic stall configurations still remain 
a subject that is hardly covered in detail; it is indispensable to at 
least insure that the transposition of the results obtained in two- 
dimensional flow is sufficiently accurate for resolving practical problems 
facing manufacturers . 

In this paper, only some aspects of dynamic stall could be touched upon 
for reasons of clarity of presention but it is certain that other interesting 
works could be mentioned. 

Moreover, papers that will be presented during this congress on unsteady 
aerodynamics will bring new results and new elements for consideration 
which will undoubtedly result in furthering-the-understanding and prediction 
of phenomena associated with dynamic'-, stall configurations. 
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